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ABSTRACT 

We have undertaken a spectroscopic search for ultra compact dwarf galaxies (UCDs) in the dense 
core of the dynamically evolved, massive Coma cluster as part of the HST/ACS Coma Cluster Treasury 
Survey. UCD candidates were initially chosen based on color, magnitude, degree of resolution within 
the ACS images, and the known properties of Fornax and Virgo UCDs. Follow-up spectroscopy with 
Keck/LRIS confirmed 27 candidates as members of the Coma Cluster, a success rate > 60% for targeted 
objects brighter than M^; = —12. Another 14 candidates may also prove to be Coma members, but low 
signal-to-noise spectra prevent definitive conclusions. An investigation of the properties and distribution 
of the Coma UCDs finds these objects to be very similar to UCDs discovered in other environments. The 
Coma UCDs tend to be clustered around giant galaxies in the cluster core and have colors/metallicity 
that correlate with the host galaxy. With properties and a distribution similar to that of the Coma 
cluster globular cluster population, wc find strong support for a star cluster origin for the majority of 
the Coma UCDs. However, a few UCDs appear to have stellar population or structural properties which 
differentiate them from the old star cluster populations found in the Coma cluster, perhaps indicating 
that UCDs may form through multiple formation channels. 

Subject headings: galaxy clusters: individual (Coma) - galaxies: dwarf - globular clusters: general 



1. INTRODUCTION 

Spectroscopic redshifts provide the most reliable method 
for establishing membership in galaxy clusters, but are 
prohibitively time consuming considering the abundance 
of faint objects in cluster fields. Given this limitation, 
studies of faint cluster galaxies often use indirect means 
such as selecting probable members by color, surface 
brightness, and morphology. Such indirect methods, how- 
ever, require assumptions about the properties of cluster 
members. Consequently, these methods are inherently bi- 
ased. Entire populations of galaxies may be left out of 
cluster member samples. 

It had previously been assumed that compact high sur- 
face brightness objects were either background ellipticals 
or foreground stars. This was due in part to a well defined 
surface brightness-magnitude relation in which dwarf el- 
liptical galaxies trend toward lower surface br i ghtnesses 
at fainter magnitudes (Caldwell & Bothun"1987 Ferguson 
k, San dage 1988; Im pey et al.||1988f Binggcli fc Cameron 
199lf |Micskc ct al. 2004a) . Only a decade ago, spectro- 
scopic surveys of the Fornax cluster, including an all-object 
2dF spectroscopic survey, revealed a new cl ass of object 
termed ultra-compact dwarf galaxies (UCDs) ( Hilker et al 



1999 Drinkwater et al. 2000 



faint objects are offset from the magnitude-surface bright 



|2000l[PhiIIipps et al.|200ip . These 
from the magnitude-surface bright- 



ness relation of dE galaxies having high surface bright- 
nesses and very small sizes. 

This discovery spurred dedicated searches for UCDs in 
other nearby groups and clusters. A large population of 
about 60 UCDs has since been established to exist in For- 



et al.[2008f|Gregg et al.)2009[ ). A similarly large popu 
tl- 



Firth 



nax (IC oudfrooij et al. 200ll Drinkwater et al.]|2005 
has recently been identified in the Hydra cluster (Misgeld 



ation 



et all2bll), while somewhat smaller populations have been 
d i 



found in the core of the Virgo (~ 25 UCDs ) (Ha§egan et al.[ 
2005| IJones et al.||2006[ [Firth et al.||2008[ ) and Centaurus 



clusters fMieske et al 



2007a[ |2009P^~^3ingle UCD wa s 



confirmed in the Dorado group (lEvstigneeva et al.|2007a[) 



Candidate UCDs have also been identified in the Coma, 



(|Adami et al. 


2009 


Madrid et al.||2010 


Wehner & Harris 


2D07; 
2008 


Mieske et al.| 


2004bl IBlakeslee & Barber DeCraatt 


Mieske et al.||2007b 


). Only a handful of UCDs have 



et al. 2009 



panions to i solated field galaxies or in poor groups ( Hau 
Da Rocha et al.| 2011 Norris fc Kannappan 



2011 1. None have been found in the poor Local or nearby 
MSTGroup ( [Chiboucas eTaI]|2009] ) environments. 

UCDs are intrinsically faint and compact objects that 
are unresolved in ground-based surveys and have magni- 
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tudes between —13.5 < My < —10.5. Effective radii for 



these objects span the range 


7 < le < 100 pc 


Mieske 


et al.|2007a 


Jones et al.|2006 


Mieske et al.|2008b 


, which 


is larger than typical globular clusters with sizes 


-2-5 



dwarf nuclei in Virgo that are compatible, [Paudel et al. 
(2010) argue that the stripping scenario is a viable option. 



pc (Larsen et al. 2001 Jordan et al. 20051 but smaller 
than dwarf ellipticals which have sizes of a few 100 par- 



sees (Drinkwater et al. 2003). At the distance of the 



Yet other researchers suggest that with a continuum of lu- 
minosities and with colors similar to metal rich globular 
clusters, UCDs may simply cons titute a bright extensio n 
of the globular cluster sequence dWehner fc Hamsl[2007| ) . 



Coma cluster, the UCD half-light radii range corresponds 
to 0.01 — 0.2 arcsec. Thus, with the Advanced Camera for 
Surveys (ACS) ~ 0.1 arcsec resolution, the larger ones are 
just resolved. The brighter UCDs also have exceptionally 
red optical colors, lying redward of the red sequence by 
up to ~ 0.2 magnitudes in Virgo and Fornax, redder than 
typical globular clusters by 0.1 magnitudes, and r edder 
than the nuclei of dE,N galaxies ( [Mieske et al.|2006[). The 
red color may be driven b y high metallicities fMieske et al. 
2006 Chilingarian et al. [2008). With these distinct prop- 
erties, iTTsTIeaFTEaFTJCDs constitute a separate class of 
object from the abundant dE and dSph galaxies found in 
clusters and groups. 

Measurements of mass-to-light (M/L) ratios have been 
made for a few of th ese obje cts and found to range from 
2-9 dMieske et al.|20 08b: Ha segan et al.|2005[[E"vstigneeva 
et al.|2007b[|Hilker et al.|2007| . The larger values could be 
evidence for the presence of dark matter but alternative ex- 
planations include: (1) inflated mass estimates produced 
by tidal heating which inc reases the veloc ity dispersion 
(IFellhauer fc Kroupa|[2006| ) or (2) bottom ([Mieske et al . 
2008aD or top heavy iMFs ([Dabringhausen et al.[[2009 
Murray[[2009| 



If UCDs originate as remnant dE nuclei via a thres 
ing mechanism they would be expected to populate cores 
of massive and dynamically evolved clusters. Compared 
to lower density environments, they would be expected in 
greater numbers and with a broader distribution within 
the massive Coma cluster. Alternatively, if UCDs are gi- 
ant globular clusters or super star clusters they should 
be associated with individual galaxies and have properties 
similar to those of the host galaxy's globular cluster and 
stellar populations. Indeed, evidence has shown UCD pop- 
ulations in Fornax and Virgo are strongly clustered, having 
smaller velocity dispersions than cluster dwarf s or even the 
central gi ant galaxy globul ar cluster systems ( Gregg et al. 
2009; Mi eske et al, 2004a; jJones et al.[[2006l [Firth eFaT 
(2008|. |Firth et al.[l2008[ ) have furthermore found that 
the rnajority of UCDs lie within the fields of the dominant 
giant galaxies in Virgo and Fornax while very few have 
been discovered to lie in intracluster regions. Therefore, 
we might expect to find a population of UCDs within the 
Coma cluster, associated with either the cluster potential 
or individual giant galaxies depending on formation mech- 
anism. 

This work is part of the larger HST/ ACS Coma Clus- 



ter Treasury Survey (Carter et al. 



2008 1 , a two-passband 



matic objects ( 


Drinkwater et al.|2000 [Fellhauer fc Kroupa 


2002; 
et al. 


Bekki et al.||2003| |Ha§egan et al.||2005t |Chilingarian 


,2008; Evstigneeva et al. 


2008 


). UCDs may simply 



A number of potential formation mechanisms have been imaging survey designed to cover 740 arcmin^ in the Coma 

Cluster to a depth of Ic ~ 26.6 mag for point sources. 
These data are being used to perform comprehensive struc- 
tural, photometric, and morphological studies of the Coma 
Cluster members. Only 28% of the originally proposed 
areal coverage was completed due to the failure of ACS, 
but it includes much of the core region. The goal of the 
Keck/Low-Resolution Imaging Spectrometer (LRIS) study 
discussed in this work was to measure redshifts and estab- 
lish membership for galaxies at the faint-end of the cluster 
luminosity function. We specifically targeted two samples 
of faint galax i es. T he first, which is discussed in [Chi 



tail of the globular cluster luminosity function or super 
star cluster end products from the mergers of massive star 
clusters themselves formed during galaxy mergers. Alter- 
natively, they may be the visible nuclei of dwarf ellipti- 
cals with exceptionally low surface brightness envelopes, 
primordial objects formed of the small scale peaks in the 
initial power law spectrum, or remnant nuclei of tidally 
stripped nucleated dwarf elliptical or late type galaxies 
which have been "threshed" during tidal encounters with 
giant galaxies and with the cluster potential. While this 
latter explanation is sometimes favored, other possibilities 
have not been ruled out. In particular, evidence exists 
showing that UCDs may be more metal rich and, in some 
cases, older than the nuclei of dE,N from which they pur- 
porte dly originated (Chilingarian et al.|2008 Mieske et al. 



2006). Because dE,N can continue to form stars in their 



cores through gas accretion while threshed nuclei, hav- 
ing lost their outer halo, suffer from strangulation and a 
cessation of star formation, it is difficult to reconcile the 
higher metallicities and, at the same time, older ages with 
a threshing mechanism. These studies instead support a 
super star cluster scenario originating from the merger of 
young massive clusters at early times. In addition, UCD 
sizes tend to be larger t han both globu lar clusters and the 
nuclei of dwarf galaxies (De Propris et al. 2005). However, 



boucas et al. ([2010[), targeted galaxies in the magnitude 
range 19 < i? < 22(-16 < Mr < -13) having low surface 
brightness and membership previously estimated through 
indirect means. The second sample, which we discuss here, 
targeted candidate UCDs with i? < 24 {Mr < -11). Can- 
didates were chosen in the core region of the Coma Cluster 
where HST/ACS data had already been obtained. 

In Section [2] we describe the observations and data re- 
duction procedure, in Section [3[ we detail the properties 
and distribution of confirmed cluster member UCDs, and 
in Section [4[ we discuss these results and potential origins. 
A summary is presented in Section [5] Throughout this 
paper, we assume a distance to the Coma Cluster of 100 



Mpc and a distance modulus = 35 (see Carter et al. 
2008| Table 1). 



finding age and metallicity measurements for UCDs and 



2. OBSERVATIONS 

2.1. Target Selection 

When this project was begun, the ACS data had just 
been taken and images were not yet calibrated. Therefore, 
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to choose a sample of potent ial UCDs, we first identified 
point sources in the catalog of Adami et al. ( |2006 ) based on 
their large ground-based Coma Cluster CFHT/CFH12K 
survey. We then imposed the following criteria: 

1. i? < 24 (Vega mag). 

2. 0.45 <{B~V) < 1.1 (3" apertures). 

3. Higher priority to sources with 0.15 < i?— / < 0.6. 

4. Located in the core region with processed ACS im- 
ages available before 2007 Feb 11. 

The B — V range includes the Fornax/ Virgo UCDs at the 
blue end and cEs at the red end. In part because of large 
errors in the Adami catalog at these faint magnitudes (the 
typical error in B — V at i? = 23.5 is 0.4 mag), we refrained 
from making our color cut too tight. 

After the initial magnitude and color selection which 
generated a list of 165 good candidates, a weight was given 
to each object based on the deviation of the ACS image ob- 
ject profile from the ACS PSF. Higher weights were given 
to those objects with any sign of having a FWHM broader 
than stellar. These weights were fed into the LRIS Au- 
toslit3 mask-making software to be used in cases of slit 
conflicts. For astrometry accurate to ~ 0.1 arcsec, we 
transformed the world coordinate system of the ACS im- 
ages to the SDSS system. Target coordinates were then 
measured from the ACS images. In an initial 2008 observ- 
ing run, a total of 47 UCD candidates were chosen by the 
software to populate 4 different masks. A second sample 
of 93 lower surface brightness galaxies were observed con- 
currently with these 4 masks with the goal of establishing 
cluster membership for faint dwarf galaxies. The results 
of this separate study are presented in [Chiboucas et al.| 
(2010) and Trentham et al. (in prep). 

The 2008 multi-object spectroscopy observations con- 
firmed 19 Coma cluster member UCDs. These original 
candidates were chosen based on the properties of known 
UCDs in the Virgo and Fornax clusters. Because the full 
range of properties for these poorly understood objects 
is not well known, we expanded our search during a 2009 
campaign to more fully explore the boundaries of the UCD 
parameter space. Targets were chosen based on a broader 
set of criteria, allowing targets with a wider range in color 
[B — V < 1.2) and with a greater range of resolution, 
including more unresolved objects. Where ACS imaging 
within the core region did not exist, we chose targets based 
strictl y on the ground-based photometry from Adami et al.| 
(2006). Greater weight was given to brighter targets, in 
particular to explore a magnitude gap (21.5 < R < 20.2) 
discovered to exist between compact elliptical galaxies and 
UCDs. 72 candidates were chosen for 2 masks, from which 
a further 8 clusters members were identified. 

2.2. Observations and Data Reduction 



Observations and data reduction are described in IChi-l 
boucas et al. (20101, but we summarize here. Because we 
are observing faint targets, down to i? < 24, along with 
low surface brightness galaxies, the spectroscopic design 
was intended to maximize light throughput at the expense 
of resolution. As we intend only to measure redshifts in 
order to distinguish between stars. Coma cluster members 



at z = 0.023, and background objects, and the Coma Clus- 
ter has a velocity dispersion of ^ 1000 km/s, a resolution 
of 200 km/s is adequate. The spectral range was chosen to 
include the 4000A Balmer break and Ca H and K lines at 
the blue end (~ 4060A at the redshift of Coma) as these 
are some of the strongest spectral features in quiescent, 
low redshift, faint galaxies. 

We use the Keck LRIS low resolution imaging spectrom- 
eter in multi-object spectroscopy mode which has very 
high sensitivity in the blue. A dichroic was used to split 
light at 5600A between red and blue chips. On the blue- 
side, we used the 400 line mm^^ grism blazed at 3400A 
providing a dispersion of 1.07A pix~^ and wavelength cov- 
erage of 4384A. With 1.2 arcsec slitlets, we achieved a res- 
olution of 7.8A FWHM. On the red-side we chose the 400 
line mm~^ grating blazed at 8500 A with wavelength cover- 
age 3950A and 1.92A pix^^ dispersion. The red-side data 
are used primarily to identify high redshift objects with 
emission line spectra. 

Four masks were observed over 2 nights, 2-3 April 2008. 
Another 2 masks were observed 28-30 March 2009 in poor 
conditions. Each mask covers a field of view of 5' x 8' with 
an average 35 slitlets. Total integrations for each mask 
during the first run were 8x 1500s with the exception of the 
blue side of one mask for which we obtained 6 x 1500s. Dur- 
ing the second run we were plagued with thick clouds and 
very poor seeing. Total integrations were ~ 8 x 1500s but 
effective exposure times are much less. Signal-to-noise A^^ 
(measured around 5000A) for the UCD candidate spectra 
range from ~ 20 to less than 1. Secure redshifts were 
measured from spectra with S/N > 4.0, while less secure 
measurements came from 2.5 < S/N < 5 spectra. Table [l] 
provides a summary of the observations. In Figure [T] we 
overlay the locations of the completed ACS fields and the 6 
LRIS masks on an image of the central region of the Coma 
cluster. The original and expanded set of candidate UCDs 
are shown as green and brown/cyan points, respectively. 

Data were reduced using the standard procedures in 
IRAF. Images were overscan corrected and corrected for 
different gains. Halogen flats for each mask were com- 
bined and a normalized flat image was generated with 
AFFLATTEN. Following division by this flat, the 8 in- 
dividual exposures for each mask were median combined 
using sigma clipping. Slit spectra were rectified by tracing 
the slit gaps for each mask and fitting these with 4th or- 
der Legendre polynomials. GEOMAP was run to compute 
the two-dimensional surface for the full set of slit gaps in 
a mask and, using that transformation, GEOTRAN was 
then executed to generate images with straightened slits. 
Arc spectra were rectified for each slitlet in the same man- 
ner. 

The usual IRAF tasks IDENTIFY, FITCOORD, and 
TRANSFORM were used to wavelength calibrate the ob- 
ject spectra from arc lamp spectra. Since arc spectra were 
taken only once per night, sky lines in each object spec- 
tra were used to correct for offsets from the lamp wave- 
length calibrations. Unfortunately, the prominent 5577A 
sky line fell on the edge of our spectra, or depending on 
the location in the mask, off the observed blue-side spec- 
tra altogether. Due to uncertainties in the applied shift, 
we therefore expect systematic errors of up to 100 km/s in 
the radial velocity measurements from our 2008 run tar- 
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gets. Observations taken in 2009 were bracketed with arc 
lamp exposures and are not expected to suffer from these 
systematic offsets. 

To extract one-dimensional spectra, we used APALL 
to identify the spectrum center, width, and sky regions. 
RVSAO/XCSAO was used, along with absorption and 
emission line template spectra, to measure redshifts. A 
flux standard was observed during a later run with the 
same configuration, used for relative fiux calibration. 

3. RESULTS 



Tables 2][3 list spectroscopic redshift measurements for 
the UCD candidate samples. We assume that objects with 
radial velocities between 4000 < Vr < 10000 km s~^, 
within 3(7 of the cluster mean, are cluster members. In 
total, we find 27 compact sources with redshift measure- 
ments consistent with Coma cluster measurement. An- 
other 14 also have measured redshifts that would place 
them in the Coma cluster, but these come from very low 
S/N (< 5) spectra and we consider these measurements to 
be highly uncertain. Out of the first set of 47 targeted can- 
didates, 19 proved to be members, 6 were highly uncertain 
members, 4 turned out to be background galaxies, and 6 
proved to be foreground stars. The remaining 12 targets 
had spectra with too low S/N to even attempt redshift 
measurements. Spectra for these 19 confirmed UCDs are 
shown in Figures [2] - |4j Thumbnail images of these con- 
firmed UCDs are provided in Fig |6] The second round 
of observations based on our expanded candidate sample 
turned up 7 background galaxies and 6 stars along with 8 
members. Spectra for these, taken in poor conditions, are 
31 and thumbnails are presented in Fig- 
Tl we display the spatial distribution of 
all confirmed and questionable UCDs along with the full 
set of candidates. Spectroscopically determined stars and 
background galaxies are also highlighted. 

In Figure [8] we show color-magnitude diagrams for all 
objects in our field. Boxes delineate the original candidate 
color selection criteria. Candidates, confirmed members, 
background galaxies, and stars are denoted by different 
symbols. dEs confirmed as members from LRIS spectra 
taken concurrently with the UCD observations are also 
shown. The dE red sequence is obvious in B — V. UCDs 
meanwhile exhibit a very large spread in colors ranging 
from the red sequence on the blue side to 0.4 magnitudes 
redder than the red sequence in B — V. One discrepant 
point with very red colors is likely affected by photometric 
errors. Two cE galaxies observed during our initial LRIS 
campaign, along with several others taken from |Price et al.| 
(2009), also lie redward of the red sequence by about 0.2 



shown in Figure 
urc [7j In Figure 



magnitudes. 

We find that, in good conditions, S/N starts to become 
too low to measure redshifts for these high surface bright- 
ness objects at i? '-^ 23.3, and we fail to measure secure 
redshifts altogether at i? = 23.5. Therefore, brighter than 
R = 23.3 where we can measure redshifts for nearly 100% 
of our sources, we find that 66% (18/27) of the original 
targeted candidates are bona fide Coma cluster members. 
This is an exceptionally high success rate given that the 
UCD candidates were selected primarily based only on 
fairly loose color criteria. Those targeted, however, were 
weighted toward objects with profiles slightly broader than 



pure PSFs and we believe that this boosted our success 
rate. 

In Figure |9j we compare the FWHM and SExtrac- 
tor classification of the UCDs as measured in the ACS 
F814W images to that of confirmed stars and galaxies. 
The ACS meas urements were perform ed with SExtractor 
as described in Hammer et al. (2010). It can be seen that 
the majority of confirmed UCDs have a FWHM at a given 
magnitude that lies between the stellar minimum and the 
larger background galaxies. One may expect that many of 
the candidates with FWHM in the range of the confirmed 
background galaxies will also prove to be background ob- 
jects, while the remainder of these points with sizes larger 
than that of the stellar minimum may turn out to be fur- 
ther UCD cluster members. These partially resolved ob- 
jects are classified by SExtractor in our images as falling 
between stars (1) and galaxies (0). Candidates in this 
range would also be considered likely members, along with 
objects classified as stars but which have slightly broader 
profiles. However, with these low number statistics, we 
cannot rule out the possibility that there are large num- 
bers of unresolved UCDs. 

3.1. Properties 

3.1.1. Photometric Properties 

In Figure [T0| we show optical color-magnitude diagrams 
for confirmed members of the Coma Cluster (UCDs and 
dEs). In all CcLS6S, cL red sequence of normal dE galax- 
ies, exhibiting some scatter, is apparent. The confirmed 
UCDs, which lie at the faint end of this sequence, have 
a much larger spread in colors of up to 0.4 magnitudes in 
both g~I and B — V. Based on the g — I ACS photometry, 
UCD colors range from just redward of the red sequence 
to nearly 0.4 mag redder. Also sho wn in these p l ots ar e 
7 confirmed Coma cluster cEs from iPrice et al.l p009k. 
These objects at brighter magnitudes also lie redward of 
the red sequence by about 0.2 magnitudes. Although there 
is a ~ 1.5 magnitude gap between the cE and red UCD 
populations, these objects do share similarly red colors 
and could perhaps form a single sequence of extremely red 
cluster objects. The colors may indicate that these objects 
have similar stellar populations and perhaps evolutionary 
histories. 

We plot UCDs and other object types in magnitude- 
surface brightness space in Figure [TT] Central surface 
brightness vs. i?— band total magnitu de comes from 
ground-based photometry ( Adami et al.| 2006 ) for all ob- 
jects in our survey region while (/i)e vs. F814W arc mea- 
sured from our ACS data. LRIS confirmed normal dwarf 
galaxies are found to lie along the well known magnitude- 
surface brightness relation for dE galaxies (I Caldwell & 



BothunI [19871 |Ferguson fc Sandage| [19881 BinggclP^ 



Cameron 1991 IBoselli et al. 2008). At much higher sur- 



face brightness, a second almost parallel sequence of UCDs 
is found, bounded at high surface brightness by the see- 
ing disk in the ground-based data where the UCDs are 
merged with the stellar sequence. The ACS data, on the 
other hand, do a much better job at separating stars and 
UCDs, although some overlap is still present. In the re- 
gion between the dE and UCD sequences are confirmed 



and probable background galaxies (see Chiboucas et al 
( )2010J ). Brightward of the UCD sequence, the set of 7 cE 
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galaxies also lie at very high surface brightness, well sep- 
arated from the normal galaxies. Two of these objects for 
which we obtained LRIS spectra are filled in. Although the 
UCDs and cEs are separated by a gap of 1.5 magnitudes, 
their surface brightness characteristics also suggest that 
these two object types may follow a continuous sequence 
and form a single class of object. 

3.1.2. Structural Properties 

Because the UCDs are marginally resolved in the ACS 
data, we are able to measure the sizes of many of these 
objects. Other studies have shown that it is possible to 
measure globular cluster effective radii if the stellar PSF 
is accurately known and the cluster size (r^) is larger than 



about 10% of the PSF FWHM (see [Harris et al.l2010[ and 
references therein). The stellar PSF for the ACS/WFC is 
about 0.1 arcsec, so we would expect to be able to mea- 
sure sizes down to about 5pc. Recently, 



Madrid et al. 



(2010) identified candidate UCDs near NGC 4874 based 



on size measurements of stellar-like objects in these ACS 
data. From a comparison of measurements in two different 
bands, they determined that they could accurately mea- 
sure sizes down to = 9.2 pc. 

To measure the structural parameters for these UCDs, 
including physical sizes, we perform surface brightness pro- 
file fitting on the objects in our F814W ACS images. Two 
methods are used to fit these barely resolved objects. We 
first model each object with a Sersic f unction using th e 
two-dimensional fitting algorithm Galfit ( Peng et al. 2002 1 . 
Galfit fits the light profile of galaxies with axisymmetric 
analytical functions and is robust for well resolved, ex- 
tended objects. Of ten used for profile fitting of globu- 
lar clusters, Ishape (Larsen 1999) is a routine in the bao- 



lab data reduction package designed specifically for fitting 
marginally resolved sources. Both convolve a PSF with 
model profiles and both algorithms determine best fits 
through minimization. 

To run Galfit, the user must provide a PSF for convo- 
lution with a model, an uncertainty map which is critical 
for determining when the best fit solution is reached, the 
choice of profile to fit, and initial guesses for each param- 
eter to be fit. The procedure we use is similar to that ex- 
plained in detail in|Hoyos et al. (2011). Briefly, to create 
the uncertainty images, we make use of the Multidrizzle 
output inverse variance images and in addition take into 
account the Poisson noise from the sources themselves. 
To generate PSFs appropriate for Galfit input, TinyTim 
(Krist 1995) PSFs are generated for a grid of locations 
on the ACS chips for each filter. These are then added 
at the corresponding shifted locations to blank individual 
distorted FLT images. The images are then run through 
Multidrizzle with the same configuration used to produce 
the final combined science images. PSF images 3 arcsec 
in size are extracted from this drizzle combined image. 
For the UCD fits, we perform the fitting in image sections 
300 X 300 pixels in size. All objects other than the central 
UCD are masked. We fit for all parameters of the Ser- 
sic function (total magnitude, effective radius Re, index 
n, axis ratio b/a, and position angle) along with the sky 
value. We found that measured values often depended on 
initial guesses. Furthermore, the Sersic index n and Rg 
are strongly correlated. We therefore ran Galfit with ini- 



tial guesses for the Sersic index ranging from 1.4 to 6.4. 
We take the fit with the lowest value as the best fit. 
Since Galfit measurement uncertainties are in most cases 
underestimates, we take for measurement uncertainties the 
standard deviation of the measurement values from fits 
with 5 sets of initial guesses, including two cases where 
we have held the Sersic index n fixed at 2.5 and where we 
included an upper limit constraint of 7.6 for the index. 

Ishape requires 10 times subsampled PSFs so we gener- 
ate an empirical PSF based on 26 real stars in one of our 
ACS images using the routines in the IRAF daophot pack- 
age. We fit each UCD with a set of 5 profile shapes: a Mof- 
fat profile with power index 2.5, a Sersic function where the 
index n is fitted for, and King profiles with three distinct 
values for the concentration parameter (defined as rt/rc 
where rt is the tidal radius and Tc is the core radius): 15, 
30, and 100. We chose to fit these as elliptical functions. 
Output includes the FWHM along both major and minor 
axes, position angle, and for the fit. To convert FWHM 
to effective radii, we take the circularized FWHM from the 
geometric mean of the semi-major and minor axes and con- 
vert to Re using the appropriate concentration -dependent 
conversion factors described in 
vided in the user guide. 



Larsen (20011, and pro- 
We take as final measurements 



the parameters from the fit producing the best residuals 
and lowest ■ These proved consistently to be from King 
models with concentration parameters 30 and 100. During 
testing, it was discovered that PSF size and fitting radius 
affected the size measurements by amounts greater than 
the small quoted errors. Since the quoted errors are there- 
fore expected to be underestimates, which also assume a 
particular model is a good match to the true profile shape, 
we therefore use the standard deviation of measurements 
from fits with 5 different models using two different fitting 
radii each (10 and 25 pixels) as a more realistic measure of 
the uncertainty in size. Final sizes are taken to be those 
based on 25 pixel fitting radii which typically produced 
the fits with the smallest residuals. 

Size measurements are presented in Table |4] To con- 
vert effective sizes from arcsec to parsec we assume a dis- 
tance of 100 Mpc. We find sizes ranging from 5 — 125 pc 
with a median effective radius (for both Galfit and Ishape 
measurements) of 23 pc. Overall, the consistency between 
Galfit and Ishape measurements of Re is quite good. We 
show the comparison in Figure \l2\ For the largest object, 
151072, Galfit finds a size about twice as large as Ishape. 
This may be an indication that the object has an extended 
halo which affects the Sersic function fits more than the 
King profile fits. Although we only fit single component 
model profiles for these objects we note that a number 
of the objects, including some of the larger UCDs, may 
be better fit with two components. For objects best fit 
with very large n, this may be caused by a lower surface 
brightness envelope surrounding a high surface brightness 
core forcing a fit with larger wings. In particular, ob- 
ject 150000, although highly obscured in the diffuse light 
of a nearby bright galaxy, upon closer inspection appears 
to have a large low surface brightness envelope. Objects 
121666 and 195526 also may require a second component 
fit. In Figure [l3j we show the residuals from both Galfit 
and Ishape fits for 4 UCDs which may be better fit with 
two components. Object 150000 in the top panel has 19% 
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of the flux remaining in the residual when fit with a sin- 
gle Sersic function. UCD 15f072 in the bottom panel has 
residuals at the level of 5% when fit with a King profile. 
Figure [M] displays residuals for 4 cases where the UCDs 
are well fit with single component models by both Galfit 
and Ishape. 

Three of the objects less well fit by single component 
profiles turn out to be the three brightest UCDs in the 
sample. If the poor fits are due to an extra envelope 
of material surrounding the UCD core, we could be see- 
ing remnant tidal debris from the stripping of nucleated, 
more massive galaxies. Extended halos have also been 
foun d around UCDs in Virgo and Fornax from HST imag- 



the Milky Way 


and other nearby galaxies (McLaughlin & 


van der Marel 


2005 




McLaughlin et al. 


2008). They may 



not necessarily evidence for the threshing mechanism. In 
fact, the lack of such signatures around most of the fainter 
UCDs in our sample may simply be due to the lower S/N 
and the small sizes of these objects. If the brightness and 
extent of the envelope is correlated with UCD size and 
brightness, extended halos around fainter objects could 
remain undetectable in our imaging. Object 150000, on 
the other hand, has an intermediate core size and bright- 
ness but more extended low surface brightness envelope, 
distinguishing it from the rest of the objects in the UCD 
sample. 

The Sersic index n measurements listed in Table |4] range 
from 1.7 - 7.6 with the majority having n > 4. Since typ- 
ically only giant elliptical galaxies are best fit with such 
large n, and since Re and n are known to be correlated, 
we compare those results to fits with forced n = 2, more 
typica l of globular clusters in the Milky Way fMcLaughlin 



et al.||2008 ). For most UCDs the resultant Re was little 
changed. Excluding 4 cases, the average difference in Rg 
between the best fit listed in the table and a fit with forced 
n = 2 is only 0.1 pc. This suggests that for such small ob- 
jects, the shape parameter has little effect on the fit, and 
cannot be well constrained. For 4 objects, the measured 
size with n = 2 was significantly smaller (163400, 151072, 
150000, and 1043225 with sizes 41.1 pc, 67.6 pc, 22.7 pc, 
and 11.9 pc, respectively), more in line with the Ishape 
size measurement. In at least 2 of these cases, the larger 
Re (and Sersic index) from the unconstrained fit is likely 
due to the presence of an outer envelope. 

3.1.3. Ages and Metallicities 

From our spectra, we measure absorption line index 
strengths in the Lick/IDS system and compare to SSP 
models to derive luminosity weighted metallicity and age 
estimates for our confirmed UC Ds. We have chosen to use 
the models of 'Schi avon ( 2007[) and the p ublicly available 
EZ-Ages code ^Graves fc Schiavon 20081) for deriving the 
ages and abundances. An advantage of EZ-Ages is that 
a separate code within this package measures line indices 
from our spectra. These models also include the effects of 
non-solar abundance patterns. 

Line index strengths from our LRIS blue-side spectra 
are measured after first degrading the spectra to the Lick 
resolution. These line strengths are fed into the EZ-Ages 
code which first determines an initial estimate for the age 



and metallicity from the II/3 and Fe indices using grid in- 
version, and then calculates the [Mg/Fe] ratio to obtain an 
alpha element abundance measurement. Age and metallic- 
ity are re-derived from H/3 vs. Mgb, where the Mgb index 
is highly sensitive to alpha element abundance. If signifi- 
cantly different from the first age and metallicity measure- 
ments, the [Mg/Fe] ratio is increased incrementally and it- 
erations proceed until a user supplied tolerance is reached. 
The process is then repeated for other Lick indices. The 
models span a wide range of age and metallicities, but will 
fail to produce derived quantities in cases where measured 
indices fall outside of the model range. 

As initial input for the fitting process, we use a Salpeter 
IMF and solar scaled isochrones. Velocity dispersions are 
also required but due to the low resolution of our spec- 
tra, we have not attempted to measure these for our 
UCDs. However, velocity dispersions have been measured 



JIIa§e 


gan et 


al.ii2005 Mieske et al.| 


2008b 


Evstigneeva 


let al. 


2007b; Hilker et al. 2007p and 


found to range be- 



twecn 9 < a < 42 km/s. We therefore simply assume a 
small ay of 20 km/s. We confirm that including ve locity 
dispersions up to 65 km/s (the mean value found by Price 



et al. ( 2009 ) for 7 cE galaxies and presumably a very high 
upper limit for our much smaller UCDs) affects derived 
values insignificantly within the uncertainties provided by 
EZ-Ages. As we did not obtain any Lick standards we 
cannot correct for any systematic offsets. Although the 
(Fe) index (an average of Fe 5270A and Fe 5335A mea- 
surements) is often used to estimate metallicity, the Fe 
5335A line is near the edge of our blue-side spectra or 
missing altogether. We therefore make use only of the Fe 
5270A line index. The strong Mgb line is used for the Mg 
measurements. 

We initially ran EZ-Ages on our UCDs individually, but 
as most have S/N < 15, we found results unreliable, hav- 
ing exceptionally large uncertainties for the derived ages 
and metallicities. We therefore stack multiple UCD spec- 
tra. For this to be useful, the objects must have simi- 
lar stellar populations and spectral properties. Otherwise, 
this would only produce a spectrum with a random mix 
of line strengths. Therefore, we stack spectra of similar 
type objects based on physical properties such as color 
and location in cluster, and on properties of the individ- 
ual line index measurements. Different sets of UCD spec- 
tra are combined. We have stacked the spectra from the 
5 brightest UCDs, 6 red (F - J > 1.05) UCDs, 9 blue 
(V - / < 1.05) UCDs, 5 red UCDs around NGC 4874, 
8 with weak II/3(< 2.5), and 8 with strong II/3(> 2.5). 
Several UCDs do not have V — I measurements and are 
not included in the color selected samples, and we do not 
include any UCDs from our second observing run as the 
S/N was particularly low. We plot II/3 and Fe 5270A line 
index measurements from the composite spectra in Figure 
1151 We also include line index measurements for 3 individ- 
ual UCDs with high S /N spectra which suggest population 
ages or abundances that are different from the majority of 
the UCDs in our sample. The uncertainties for these are 
large, however, with error bars which span much of the 
grid. For clarity, we therefore do not include error bars for 
these individual objects. Overlaid in the figure are grids 
for a/Fe = and 0.3 with lines of constant age and [Fe/H]. 
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Spe ctra for non-UCDs come from |Chiboucas et aL] 
(2010). Included in this figure are measurements for ellip- 
tical, dE, and dE,N galaxies. To produce these, we have 
stacked spectra for 2 small elliptical galaxies, 10 low sur- 
face brightness dE galaxies, 4 dE,N with prominent nuclei 
and 24 dE,N with small/faint nuclear star clusters. For the 
4 objects with prominent nuclei, we re-extracted the two 
dimensional spectra with smaller apertures meant to in- 
clude light primarily from the nucleus although some con- 
tamination from the spheroidal component of these dE,N 
may still be present. We have also measured Fe 5270 A and 
H/3 for one of t wo cEs for which we hav e obtained spectra 
(object 91543 in Chiboucas et al. (20101). This object was 
also observed with Hectospec o n the MMT and lin e index 

poool). 



Both 



measurements are presented in Price et al. 
measurements are included in this plot. 

Table [5] presents line index measurements and derived 
ages and metallicities ([Fe/H]), along with [Mg/Fe], an in- 
dicator of alpha element abundance. For comparison, we 
include in this table [Fe/H] estimates from the V — I col- 
ors using an e mpirical relation based on Galactic globular 
clusters from Barmby et al. (20001. For the composite 
spectra cases, we use an average V — I color. 

Given the large uncertainties, it is clear that our age 
and metallicity measurements should not be taken at face 
value. Furthermore, we find that the addition of single 
UCDs in composite spectra in some cases produce a large 
change in measured values. However, we can draw some 
conclusions from this exercise. The most obvious is that 
the UCDs are not a homogeneous population. Rather, 
they exhibit a spread in age and metallicity significant at 
the > 2.5(7 level. This was already expected based on the 
large color spread we find for these UCDs. Line index 
measurements confirm this range is due to real differences 
in stellar populations. Redder UCDs are found to be both 
older and more metal rich than the bluest UCDs. Blue 
UCDs have on average intermediate ages and low metal- 
licity. The majority have sub-solar metallicity, although 
a couple individual cases show evidence for more metal 
rich stellar populations. Red UCDs around NGC 4874 are 
found to be no different from the overall set of red UCDs. 

With the exception of a subset of red UCDs, most UCDs 
appear to have high [Mg/Fe] abundances, suggesting pos- 
sible super-solar alpha element abundances. This would 
indicate that the stellar populations in UCDs formed in 
quick bursts of star formation. Because of the large un- 
certainties for this derived quantity, we caution the reader 
against drawing any firm conclusions. However, we note 
with interest that the derived values for most subsets of 
UCDs are similar to those of globular clusters which typi- 
cally have super-so lar alpha element abundances. See e.g. 
Puzia et al.| ( [20051 ) who find a mean [a/Fe] = 0.47 ± 0.06 
with a dispersion of 0.26 dex for globular clusters in early 
type galaxies. The set of strong H/3 UCDs are found to 
be very young, metal poor, and highly alpha element en- 
riched. The relation between star formation time scale 



and alpha element abundance in Thomas et al. (2005) 
would suggest a very rapid burst of star formation with 
a timescale of only ~ 0.25 Myr. The extreme values for 
this set of objects is likely strongly influenced by the in- 
clusion of object ID 150000. This object appears to be 
surrounded by a very faint extended envelope and may 



not be the same class of object as these other UCDs. 

Comparing dE to dE,N derived ages and metallicities, 
we do not find any differences within the uncertainties. 
Red UCDs are generally older than the dE and dE,N galax- 
ies, significant at greater than 2a. There is a slight hint 
that the UCDs are more metal rich as well, but this is 
insignificant within the errors. Blue UCDs are similar to 
the dE and dE,N within the errors, having slightly lower 
metallicities and similar ages. The measured indices we 
find for the one cE are quite different from the UCDs, 
having super-solar m etallicity and a you ng to intermediate 



age. Other cEs from 

as low as those found for the red UCDs 



3.2. Distribution 



Price et al. ( |2009 ) have metallicities 



Looking at the spatial distribution displayed in Figure 
[TJ we find that all confirmed UCDs lie toward the cen- 
tral core region of the cluster along a band in declina- 
tion a.t S ^ 27.97. We do not find any members lying far 
from this, although it can be seen that very few UCD 
candidates were targeted away from this band. Three 
crosses (denoting background galaxies or foreground stars) 
lie north of 28.06 and another three south of 27.94. With 
these small number statistics it is hard to say too much 
about the true spatial distribution of the full UCD pop- 
ulation. However, if we look at the original set of good 
candidates (from which we achieved a 66% confirmation 
success rate) in conjunction with the confirmed members, 
we find a pronounced linear structure running E-W from 
NGC 4874, past the other central giant eUiptical NGC 
4889, and continuing toward IC 4051. For a comparison, 
we take a sample of 10,000 points distributed randomly 
over the fields which were available during the initial se- 
lection of candidates, assuming a uniform distribution. A 
2-dimensional Kolmogorov-Smirnov (KS) test finds that 
the good candidates have a distribution different from a 
uniformly distributed sample at a 99.95% confidence while 
a 1-dimensional KS test finds that the candidate popula- 
tion differs in declination at 99.997% confidence. In con- 
trast, the candidates have a different distribution in RA 
than the random sample at only 67.18% confidence. As it 
is hard to imagine how the color selection criteria could 
have biased the candidate list spatially, we suspect this is 
a real structure in the Coma cluster core. It is possible, 
however, that this apparent linear structure is produced 
primarily by UCDs associated with giant galaxies in the 
cluster core. 

We also see what appears to be a concentration of con- 
firmed UCDs around the cD galaxy NGC 4874. A ma- 
jority 7 of 9 objects around NGC 4874 have red colors 
(V — I > 1.05), while UCDs east of this grouping have 
a greater percentage of blue V — I colors (at least 10 of 
16 objects). We also find that aU UCDs near NGC 4874 
have Vr > 6800 km s~^. Objects eastward of this display 
a larger range in measured radial velocities. 

Turning to the radial velocity distribution, we note that 
line-of-sight velocity dispersions for UCD populations in 
other clusters have been found to be smaller than for other 
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populations of 979 ± 30 and 1096 ± 45 km s~ , respec- 
tively. The Coma cluster has a mean radial velocity of 
6925 km s~^. For the faint, low surface brightness sample 
of 51 dwarf cluster members observed with LR IS, we find 
(vr-) = 6970± 178 with cr^, = 1269±126 km s^^ (jChiboucas 



|et al.|[2010| , slightly higher than previous measurements 
For the UCD population, we find a mean radial velocity of 
(vr) = 6887±207 and a velocity dispersion (Ty = 1072±146 
km s^^, comparable to what had been found for giant and 
brighter dwarf galaxies, and smaller (although by < 2a) 
from what we find for dE and dE,N galaxies in the same 
region. Histograms showing the radial velocity distribu- 
tion of the UCDs along with that of the normal dwarfs 
observed in the same masks are presented in Figure 16 Pe- 
culiar velocities of prominent galaxies are indicated. The 
left histogram shows that most UCDs have velocities sim- 
ilar to, and centered on, the cluster mean, with just a few 
outliers. 

The isolation in color, space, and velocity noted in the 
spatial distribution is also evident in the radial velocity 
histograms (Figure 16 ) where there appears to be a distinc- 
tion between the UUD groups split either by color or RA 
with hints of associations with specific major galaxies. For 
example, 12 UCDs with RA < 195 have (vr) = 7296 ± 160 
km with = 558 ± 119 as compared to 15 ob- 
jects at RA > 195 with (vr) = 6559 ± 325 km s'^ and 
(7v = 1257 ± 238, a difference of about 3a in the velocity 
dispersion. Table [6] lists mean radial velocities and disper- 
sions for different UCD subsamples. 

Putting these distribution trends together in Figure [iTj 
we show UCD projected distance from 3 prominent Coma 
giants vs. radial velocity and see evidence that many of 
the UCDs are associated with giant ellipticals. In particu- 
lar 9 UCDs in velocity and spatial proximity to NGC 4874 
have (vr) = 7257 ± 87 with cr„ = 254 ± 60 km s"^ This 
is very similar to the radial velocity, Vr = 7220 km s~^, 
of NGC 4874 and over 3a from the cluster mean, indicat- 
ing these objects are more likely to be associated with the 
cD galaxy NGC 4874 than the general cluster potential. 
A stud y of Coma cluster globular clusters by [Peng et al.| 
( [20TT| ) finds that the GC population of NGC 4874 extends 
out to ~ 130 kpc before the intracluster GC population 
starts to dominate. This corresponds to about 4.5 arcmin. 

We calculate whether the 9 UCDs near NGC 4874 are 
likely to be bound to the cD galaxy. The mass of NGC 
4874 has b een estimated at 1.4 x 10^^ Mq from X-ray ob- 
servations ( |Vikhlinin et al.||l994[ ). We take 1.0 x 10^^ Mq 
as a lower limit. The escape velocity for NGC 4874 is then 
given by 

Vesc = y(2GM/r,ep). (1) 
For the UCD at the largest separation from NGC 4874 
(121666, 3.7 arcmin), we find v^sc — 863 km s~^ at a pro- 
jected distance of 108 kpc. The difference in radial veloci- 
ties for the 2 objects is 369 km s~^, lower than the escape 
velocity a nd therefore consist ent with being a bound satel 
hte. As 



Firth et al. (2007), we confirm that this object 



also lies within the tidal radius of NGC 4874 as imposed 
by th e othe r central giant galaxy NGC 4889. |Vikhlinin| 
etaTl ((19941) find similar masses for NGC 4889 and NGC 



4874. 



and M are the masses of NGC 4874 and NGC 4889, re- 
spectively, and D is the separation, we find a tidal radius 
of ~ 5 arcmin. All 9 UCDs are within this projected dis- 
tance, and with a velocity dispersion of only 254 km s^^, 
we expect these UCDs to be bound to NGC 4874. 

Although clustering in both velocity and spatially 
around NGC 4889 is not as clear, the 5 UCDs in closest 
proximity, within 4.5 arcmin of NGC 4889 at Vr = 6495 km 
s-\ have (vr) = 6526 ± 148 km s^^ with ay = 332 ± 105. 
This is also nearly 3a from the cluster mean. The lack of 
confirmed UCDs within 2 arcmin of this giant elliptical is 
strictly due to the fact that the ACS field centered on this 
galaxy was never observed. 

Three very low velocity UCDs which are spatially coin- 
cident with IC 4051 have (vr) = 4706 ± 66 km s'^ with 
ay = 114 ± 47, similar to the very low peculiar velocity 
of 4779 km of this giant. This galaxy, which is not a 
central dominant giant, is remarkable for having one of the 
highest globular cl uster specific frequencies (Sat = 1 2.7) in 
the Coma cluster |Marm-Franch fc Aparicio 2002[ ). Two 
of the UCDs have a fairly large projected separation from 
this galaxy. It is possible that while these UCDs may not 
be bound to this giant, they may be part of a cluster sub- 
structure which includes IC 4051. 

In addition, two objects with particularly large radial 
velocities are found within 4 arcmin of IC 3998 having 
a similarly large radial velocity. We also identify several 
UCDs which are coincident with the galaxies IC 4042 and 
IC 4041 toward the eastern side of the cluster. As the sys- 
temic velocities for these two galaxies are similar to that 
of the cluster mean, it is possible that these UCDs simply 
belong to the general cluster potential. Searching for all 
possible associations of UCDs with giant galaxies, we find 
nearly all UCDs have potential host giant galaxies (Figure 



17) 



rom 



rt = {m/3Mf'^D (2) 
( Binney fc Tremaiiie||1987 ) where rt is the tidal radius, m 



In the left panel of Figure [T8j we show the cumulative 
distribution of confirmed UCDs as a function of distance 
from the nearest of one of the 3 brightest galaxies in the 
core region (NGC 4874, NGC 4889, and IC 4051). For 
comparison we determine the expectation for a uniform 
distribution. This is done by randomly distributing 25,000 
points over the LRIS footprint and normalizing the resul- 
tant minimum distance cumulative distribution. We find 
that the UCDs are more concentrated around the giant 
galaxies than a spatially uniform distribution would pre- 
dict. The right panel is similar but includes velocity in- 
formation. In this case, the cumulative velocity difference 
from the nearest of one of the 3 giants is shown. Velocities 
are attached to the random sample of 25,000 points by 
assuming a Gaussian distribution for the velocities having 
(w,.) — 6925 km/s and a = 1000 km/s. Again we find that 
the UCDs are more strongly clustered around the giant 
galaxies. A clear 2/3 of the confirmed UCDs are associ- 
ated with one of 3 central giants. 

With the exception of the 2 central cD/giant ellipticals, 
other UCDs may not be bound to their respective neigh- 
boring galaxies. However, if this is the case, the velocity 
and location groupings we find suggest that either these 
are previously bound clusters in the process of dissolution 
to the cluster potential, or that they belong to substruc- 
tures present in the Coma cluster. 



UCDs in Coma 



9 



4. DISCUSSION 

From an initial observing run we identified 19 members 
selected primarily based on color along with some informa- 
tion on profile size. For targets with R < 23.3, our com- 
pleteness limit, we had an unexpectedly high 66% success 
rate for these initial observations. A follow-up run with 
an expected lower confirmation rate due to intentionally 
loose selection criteria confirmed a further 8 UCDs from 
poor weather observations. These results suggest that the 
Coma cluster harbors a large population of UCDs, at least 
in the core region, although this high success rate must be 
due at least in part to the fact that we can marginally re- 
solve these objects with the superb resolving capabilities 
of HST/ACS. 

The spatial distribution of UCDs are compared to that 
of the Coma cluster globular cluster and dE,N populations 
in Fi gure 19 Candidate globular clusters (Peng et al. 
2011 1 are mapped as points in the top plot^ There is 
a clear concentration of globular clusters visible around 
NGC 4874. An overdensity is also apparent surrounding 
the missing ACS field that would have included NGC 4889. 
On the eastern side of the cluster, there is a slight over- 
density at S ~ 28.02 deg near IC4051 which lies just off 
our ACS footprint. Intracluster globular clusters not asso- 
ciated with individual galaxies are spread throughout the 
core region. After smoo thing the globular cluster distribu- 
tion, Peng et al. (2011 1 find an excess which forms a band 
through the core between NGC 4874 in the west to IC 
4051 and NGC 4908 in the east. Similarly, the UCDs con- 
firmed to date exhibit a concentration around NGC 4874 
along with an east- west flattened distribution through the 
core forming a linear structure in the same region as this 
apparent band of globular clusters. Considering candidate 
UCDs lends further support for the presence of a band- 
like structure running through the core. However, with 
the small number statistics for UCDs, it is difficult to tell 
if this is a real structure, or if the UCDs are simply follow- 
ing the distribution of more massive galaxies in the core 
region. 

Meanwhile the dE,Ns exhibit a slight excess near NGC 
4874 and NGC 4889, but otherwise appear to be spread 
uniformly through the ACS footprint, with no indication 
of either an excess or deficit in the eastern side of the core 
where the remainder of the UCDs have been confirmed. 
In Figure |20[ we show the cumulative distribution of con- 
firmed and candidate UCDs, dE,N, and GCs as a function 
of projected distance from NGC 4889. For comparison, we 
take a sample of 10,000 points distributed randomly over 
the observed core ACS footprint. The UCDs are more cen- 
trally concentrated than any other population. The dE,N 
are also more concentrated than expected for a uniform 
distribution, while the globular clusters have a shallower 
distribution until ~ 8 arcmin in distance from NGC 4889 
most likely due to the large concentration around NGC 
4874. However, taking the cumulative distributions as a 
function of distance in declination only, we find that both 
GCs and UCDs are much more concentrated toward the 
central declination, while dE,N follow the uniform distri- 
bution. This difference between the populations is greater 
when considering declination only, most likely because of 
the strong clustering of GCs and UCDs around the central 
giants, and perhaps due to the excess of compact systems 



noted along a band in declination. 

In future work, we plan to do more detailed photometric, 
structural, and distribution comparisons between Coma 
UCDs and dE,N gala xies. For now, we u se the luminos- 
ity relation found by Cote et al. (20061 between Virgo 
dE,N and their respective nuclei to compare the magni- 
tudes of UCDs with those expected for the Coma dE,N 
nuclei. We find that within the Coma footprint, there 
are only ^ 20 — 30 dE,N within the magnitude range 
15 < F814Vl^ < 20 which could host nuclei with mag- 
nitudes in the range 21 < FSIAW < 23.5. We have con- 
firmed 27 UCDs and identified over 100 other candidates. 
We expect the ratio of UCDs/N (21 < F8UW < 23.5) 
to be at least 2-3. If we assume that all UCDs are pro- 
duced from threshing of dE,N, the efficiency of stripping 
within the Coma core region must be very high. However, 
given the requirement for highly eccentric orbits, t hresh- 
ing models do not predict such high efficiency rates ( Bekki 
et al.|2003 ). In this scenario, one might also expect to find 



the number of dE,N to be depleted toward the core of 
the cluster. Instead, the dE,N are found to have a fairly 
uniform distribution with perhaps a slight excess in the 
central core. 

The confirmed UCDs display a wide range in color ex- 
tending to exceptionally red objects. One explanation for 
the extreme red colors is that a dE,N which initially lies 
along the red sequence undergoes threshing and is stripped 
of at least 4 magnitudes of material by the tidal field of 
the cluster without affecting the color of the remnant nu- 
cleus. We find evidence for an alternative explanation 
when we compare UCD and globular cluster colors. Glob- 
ular clusters exhibit a well known bimodal color distri- 



bution likely due to a metallicity bimodality (West et al. 
2004 1. This generates a large range in globular cluster 
colors. In Figure [21] we display the g — I color distribu- 
tion of the confirmed UCDs with that of the Coma cluster 



globular cluster candidates from Peng et al. (2011 ). Small 
circles represent globular clusters from visit 19, the field 
which includes NGC 4874. A histogram of globular clus- 
ters fainter than / > 24.7 contains only one obvious peak 
around g — I = 0.95. If we look only at visit 19 globular 
clusters brighter than / < 24.7, the second peak becomes 
apparent. The histograms in this figure are scaled to fit 
within the boundaries of this plot. The total number of 
globular clusters represented by this lower histogram is 
~ 500. We fit a double Gaussian to the lower histogram 
and find peaks at 0.93 an d 1.17 with a = .09 for both. 
For the full sample of GCs, Peng et al.] (|2011|) find peaks at 
- 0.9 and 1.15. For inner GCs around NtSCJ 4874 they find 
a slightly redder blue peak at 0.94. The UCDs are found 
to span the same wide range in color as the GCs. With so 
few confirmed UCDs, it is difficult to determine whether 
they also exhibit a bimodality in color. To calculate aver- 
age colors, we assume the Gaussian distributions for the 
globular clusters in order to assign a weight to each UCD 
for grouping with the red and blue distributions. We find 
weighted averages of blue and red objects oi g — I = 0.95 
and 1.15, very similar to those of the globular cluster pop- 
ulation. Since we are assuming the same distribution, we 
cannot argue that they share this bimodality. However, 
this does indicate that the UCD distribution is not incon- 
sistent with that of the globular clusters. 
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One secure conclusion that we can draw from this plot 
is that there is no discontinuity between the two popu- 
lations; the UCDs simply extend to brighter magnitudes 
while exhibiting the same spread in col or as GCs. In fact , 
due to the lack of any discontinuity, thelPeng et al.|(|2011|) 



globular cluster catalog includes many of our confirmed 
and candidate UCDs. It is worth noting that the red GC 
peak is more prominent when considering bright GCs into 
the range of UCDs in the field around NGC 4874, where 
we also find a large fraction of red UCDs. The possible 
extension of globular clusters to brighter magnitudes has 
also been observed for the cD galaxy NGC 3311 in Hydra 
where the red globular cluster sequence is seen t o extend 



upward in luminosity into the range of UCDs (Wehner 
[etaL][2008l) . 

Both NGC 4874 and IC4051 harbor large globular clus- 
ter populations, with Sjy ~ 9.0 and 12.0, respectively 



(Marm-Franch & Aparicio 



2002 Harris et al. 2009). We 



find a large population of UCDs that is almost certainly 
bound to NGC 4874 and another group of three UCDs 
likely to be associated with IC 4051. Unfortunately, this 
latter galaxy lies east of our ACS footprint. Confirma- 
tion of a large UCD population around this non-central, 
non-cD galaxy would lend strong support for a star cluster 
origin. 

We test whether the UCDs can be accommodated by 
the bright tail of the NG C 4874 globular cl uster distri- 
bution. Gaussian fits to Peng et al. (2011) GC candi- 
dates in visit 19 find Ipeak = 26.2 with a = 1.2. How- 
ever, incompleteness sets in before the peak. We there- 
fore take the globular cluster luminosity function (GCLF) 
values found rec ently for M87 in Virgo, lAB,peak = 26.9 
with cr = 1.37 ( |Peng et al.1|2009[ |2011[ ). This produces 
fewer counts than UCD candidates at the very bright end 
by > 3(7 for F814M^ < 22. Fainter than F814I^ > 22 
{Mpiiw > —13) the UCD population is fully consistent 
with being drawn from a Gaussian distribution of GCs. In 
fact, many of the UCD candidates are included in the GC 
candidate list. We also take the M87 form of the GCLF 
and randomly distribute in magnitude according to this 
Gaussian distribution the ^ 3000 globular cluster candi- 
dates found in the same ACS field as NGC 4874. From 
1000 simulations, we find that the brightest expected glob- 
ular cluster is Mj = —13.0 ±0.4. The brightest confirmed 
UCD in this field has AIpsuw = —13.2, consistent with 
the expectations. 

We compare the sizes of the Coma cluster UCDs with 
other compact objects in Figure [22j Globular clusters in 
the range — 10 < My < — 8 have nearly constant half-light 
radii ~ 2.5± 1.5 pc independent of luminosity. Faint ward 
of this, the observed sizes increase, with effective radii as 
large as 20 pc being found for clusters as faint as —5. 
Brightward of My — —10, the few very bright globular 
clusters start to increase in size with increasing luminosity. 
These objects merge seamlessly into the range of the so- 
called UCDs which display a trend of increasing size and 
decreasing surface brightness with increasing luminosity. 
The Coma UCDs follow the same luminosity-size relation 
as found for Virgo and Fornax UCDs. Compact ellipticals 
and normal dEs are found toward the bright extension of 
this sequence, although cEs would fall off toward higher 
surface brightnesses while dEs have lower surface bright- 



ness. Nuclei of early type galaxies in the Virgo cluster 
tend to have smaller siz es at a given magnitude than the 
UCDs dCote et al.|[2006t . This was also noted for nuclear 
star cl usters in low mass dwarf galaxies by [Georgiev et al.| 
([2009]), who suggest that nuclei may expand in size upon 



being freed from the strong gravitational potential within 
a galaxy, e.g. in the case of galaxy threshing. 

A transition between objects which have sizes indepen- 
dent of magnitude and those which follow a magnitude 
- size relation has been previously noted at My ^ —11. 
This transition occurs at the sa me magnitude as a b reak 



in the metallicity distribution (Mieske et al. 2006). It 
furthermore corresponds to a mass of ~ 2.5 x lO^M© 
where a break in the log cr — logM scaling relation has 
been noted between globular clusters and larger systems. 
Higher mass-to- light ratios, in the range 6 < M/Ly < 9 
have been found for objec ts above this transition magni- 
tude ( Ha§egan et al.|[2005 ). As these cannot be explained 



with canonical IMFs and baryonic matter, they suggest 
the fundamental difference between globular clusters and 
UCDs is the presence of dark matter. This supports the 
threshing model in which remnant cores of threshed dE,N 
retain some of their dark matter halo. The nuclei them- 
selves may have formed from the coalescence of gl obular 



et al. 


2004 


Lotz et al. 


2002 



This would produce com- 
plex stellar populations and a range of colors/metaUicities, 
along with larger sizes and masses consistent with what 
is found for massive globular clust ers and perhaps UCDs 
(iGeorgiev et al.||2009l ]Hilker||2006|). 



et al. 



However, another possibility for the larger sizes and ve- 
locity dispersions of the UCDs is a diffe rence in physics 
at the time of formation. Murray] ( 2009 ) shows that star 
clusters above M > IO^Mq would be optically thick to 
IR radiation at the time of formation, and that the bal- 
ance between radiation pressure and gravity sets up a size 
- mass relationship. If the break in the IMF is set by the 
Jeans mass, they argue that optically thick clusters will 
have top heavy IMFs. Top heavy IMFs could be the ori- 
gin of the high M/L ratios, since massiv e stars will age 
quicker and leave more stellar remnants ( Dabringhausen 
^ ^[2009| . 

Taken together, the properties of the Coma cluster 
UCDs appear to point toward a star cluster origin for a 
majority of these objects. The UCDs exhibit strong spa- 
tial and velocity correlations with the major galaxies in 
the core. In particular a large fraction of the UCDs re- 
side within the halos of the massive galaxies, at least in 
the case of the two cD/giant ellipticals in the cluster core. 
Other UCDs may be associated with some of the other gi- 
ant galaxies in the core region. The UCDs are furthermore 
correlated in color and metallicity with the host galaxy. 
We find a large population of red UCDs around NGC 4874 
and a bluer population around NGC 4889. Unfortunately 
due to the missing ACS field, we do not know whether a ra- 
dial gradient in color/metallicity could be present around 
NGC 4889. The range of UCD colors is identical to that 
of globular clusters and there is no evidence for disconti- 
nuity between the two populations in luminosity. The spa- 
tial distribution shows remarkable similarities to globular 
clusters with a large number of confirmed and candidate 
UCDs found around the same galaxies hosting large glob- 
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ular cluster populations. Both UCDs and globular clusters 
also show some evidence for a structure of compact objects 
running east-west through the core region. 



From the age-metallicity diagram (Figure 15), we find 
that many of the UCDs are not extremely old, and thus not 
primordial in origin. Red UCDs are more metal rich than 
blue ones and must have undergone self-enrichment or 
been formed from pre-enriched gas. Red, metal-rich glob- 
ular clusters are more often found aro und massive galaxies 

The luminosity of 



(Peng et al. 


2006 


Harris et al. 


2009 



globular clusters in a system also correlates with the lu- 
minosity of the host galaxy, with brighter ga laxies hostin g 
populations with brighter globular clusters ( Hilker|[2009 1 . 
If UCDs simply extend the globular cluster sequence to 
brighter magnitudes, it then follows that UCDs should be 
found around the brightest and most massive galaxies such 
as the central Coma cD/giant elliptical galaxies and should 
include a significant number of red objects, as we find. 

The commonality in UCD color/metallicity properties 
with location could be a consequence of the formation of 
UCDs (and GCs) in a small number of discrete star for- 
mation events, each event characterized by the metallicity 
of gas available for star formation. These star formation 
events could be induced by e.g. cataclysmic wet mergers 
at early times. Since each galaxy has a unique history 
it follows that the UCD and GC populations should vary 
from galaxy to galaxy. The affiliation by color with host 
suggests that the UCDs were born in a single event, or 
only a couple events. 

Some of the blue UCDs are associated with NGC 4889. 
Others may follo w the globular clus ter intracluster distri- 
bution noted by Peng et al. (2011). Some of these blue 
intracluster globulars may be metal poor clusters stripped 
from infalling dwarf galaxies. In this scenario, it is pos- 
sible that the largest compact stellar systems may be the 
remnant nuclei of som e of th ese disrupted dwarf galaxies. 
However, [Peng et "aL (2011) find ^ 20% of intracluster 
globular clusters are red and suggest that at least part of 
this intracluster population may be stripped GCs from the 
halos of massive galaxies. In this scenario, the UCDs could 
also simply be giant globular clusters stripped from these 
more massive galaxies. 

The findings so far cannot rule out threshing as an ori- 
gin for UCDs. The threshing model predicts that UCDs 
will be found in greater abundance and have a distribu- 
tion concentrated toward the center of th e cluster potential 
and around cluster super giant ga laxies ( Bekki et al.|2003 



Bekki 2007 Thomas et al. 2008). This is generally what 



we find, although detailed modeling would be needed to 
understand the E-W linear structure of compact objects 
if real. The differentiation of color by host is more diffi- 
cult to explain within the framework of a threshing sce- 
nario. One possible explanation is if the color variations 
we find are due to radial gradients in metallicity in which 
more metal rich dE ,N are found closer to the central giants 
( Mieske et al.|2006 ), thereby producing an excess of metal 
rich UCDs with smaller orbital radii around the central gi- 
ants. Another possibility for explaining the large fraction 
of red UCDs around NGC 4874 is if more massive (and 
by implication, more metal rich) galaxies require smaller 
pericenter radii and interaction with more massive galax- 
ies for efficient stripping. Since the sizes of dE,N nuclei 



are known to scale with the brightness of the host dwarf 
galaxy, we would then expect that larger UCDs would be 
redder and that these would be found preferentially close 
to the central giants. From Figure [23] we see only hints 
of such trends. After two outlying UCDs are excluded, 
we find that larger objects tend to be redder, but with 
a correlation coefficient of only 0.36. We also find that 
larger objects tend to be closer to NGC 4874, but this 
trend is also very weak, with a correlation coefficient of 
-0.26. No correlation is found when comparing UCD size 
with distance from NGC 4889, or with the distance from 
the closest of either NGC 4889 or NGC 4874, although 
the missing ACS panel centered on NGC 4889 may be the 
cause of this. 

Some similarities between the UCDs and cE galaxies, in 
particular in the color and surface brightness characteris- 
tics, are seen. There may be some correlation in spatial 
distribution as well, although larger numbers of cEs would 
be needed to establish this. Since cE galaxies are expected 
to form from tidal stripping of more massive galaxies, a re- 
lation between the two populations would suggest a com- 
mon threshing origin. However, we have thus far been 
unable to fill a ~ 1.5 mag luminosity gap that currently 
separates these two object types. 

We have identified a few confirmed UCDs which may 
have very faint, low surface brightness envelopes. One of 
these, object 151072, has obvious surrounding structure, 
the largest measured size of our confirmed UCDs, and a 
fairly blue color. Another object, 150000, has a hint of a 
large very low surface brightness envelope and a young de- 
rived age. Since it is expected that destructive processes 
such as threshing should play a role in rich clusters like 
Coma it might be expected that at least some fraction of 
these objects are remnants from galactic stripping. |Nor-| 
ris & Kannappan pOll ) make a statistical argument for a 



bright upper limit for star cluster formation. They suggest 
that compact objects with My > — 13 may be a mixture 
of objects formed through star cluster formation processes 
and those formed via a stripping mechanism. Compact 
objects brighter than this limit cannot be formed as star 
clusters or through the merger of star clusters simply be- 
cause there are no globular cluster systems which are abun- 
dant enough to be populated that far into the bright tail 
of the GCLF. Object 151072 is our brightest UCD with 
My — —12.9, close to this suggested limit. Two other 
very bright UCDs also display hints of surrounding struc- 
ture. Object 150000, while not particularly bright, does 
not resemble Coma cluster GCs in stellar population or 
structural properties. Thus, although our findings sug- 
gest that the majority of Coma cluster UCDs have a star 
cluster origin, several of the confirmed UCDs which have 
slightly different properties may in fact have a galactic 
origin. Indeed, a number of other recent studies have also 
found eviden c e for multiple formation channels for UCDs 
iHilkerl '2 0061 |Da Rocha et al.||2011[ [Taylor et al.||2010[ 
Chilingarian et al. 2011 Norris & Kannappan 201ip, all 



finding both cases where a star cluster origin is preferred 
and cases where stripping of a more massive galaxy is the 
more likely origin. 

The results presented here are based on UCDs located 
within the core of the Coma cluster. In future work it will 
be important to search outer regions as well, to completely 
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define the spatial and velocity distribution and range of 
UCD properties. Future work will discuss the full can- 
didate population and compare in greater detail to the 
properties of Coma cluster dE,N. 

5. SUMMARY AND CONCLUSIONS 

Using LRIS on Keck, we have spectroscopically con- 
firmed 27 UCDs within the core region of the Coma Clus- 
ter. With an initial UCD confirmation success rate greater 
than 60% for Mpsi4W < —11.7, we believe that the rich, 
dense, evolved environment of the Coma Cluster core hosts 
a large population of UCDs. However, the high suc- 
cess rate is also due in part to the fact that UCDs are 
marginally resolved in ACS imaging. We find properties of 
the UCD population consistent with what has been found 
for UCDs in other clusters, having the same range in lumi- 
nosities, sizes, colors, and metallicities and having similar 
distributions. 

The confirmed Coma cluster UCD population has mag- 
nitudes Mj > —13.5, colors 0.8 < [g — I) < 1.3, sizes 
primarily in the range 7 — 40 pc, and metallicities —1.3 < 
[Fe/H] < —0.6. They are distributed through the central 
core between NGC 4874 in the west and IC 4051 in the 
east. We find strong spatial and velocity correlations with 
the major cluster galaxies. A subset of UCDs is almost cer- 
tainly bound to cD galaxy NGC 4874. Other UCDs may 
be associated with the other central giant, NGC 4889, and 
with some of the other giants in the core region. Notably, 
the UCDs also exhibit color/metallicity correlations with 
location in the cluster. NGC 4874 hosts a large population 
of red, metal rich UCDs while NGC 4889 appears to host 
primarily blue UCDs, although a radial gradient cannot be 
ruled out. There is also a subset of blue UCDs which may 
lie in the intracluster region. The affiliations with host by 
color suggests formation in discrete star formation events 
with metallicity determined by that of the gas pool avail- 
able during such an event, e.g. during the same early time 
cataclysmic wet mergers that are believed responsible for 
producing globular cluster populations. 

We suggest that these objects could be related to glob- 
ular clusters. Not only do they share similar colors and 
lie along the extension of the bright tail of the GCLF, but 
they also have a similar distribution in the cluster. NGC 
4874 has a very high GC specific frequency and hosts a 
significant UCD population as determined by the large 
number of both confirmed and candidate UCDs surround- 
ing this galaxy. IC 4051 is another interesting galaxy lying 
east of our ACS survey region having an unusually high 
GC specific frequency and a non-central location in the 
cluster. We have identified several UCDs associated with 
this galaxy and find a slight excess nearby in our full candi- 
date list. Confirmation of a large UCD population around 



IC 4051 would provide further evidence for a GC-UCD re- 
lationship. Intracluster GCs have been found distributed 
throughout the core region with a possible excess running 
in an E-W band. We find a similar structure with our 
confirmed UCDs. Compared to dE,N, possible progenitors 
in a threshing scenario, the UCDs are more concentrated 
around massive galaxies while dE,N are dispersed through 
the core. 

Although most UCDs appear to be oversized globular 
clusters, a few of these objcicts c;xhibit more evidence for 
a threshing origin. With hints of diffuse surrounding en- 
velopes, bluer colors, younger ages, and/or exceptionally 
large sizes, these objects may in fact be remnant nuclei 
from stripping of dE,N or more massive galaxies. Since de- 
struc-tive processes are rampant in dense clusters, it would 
almost be surprising if such stripping never occurred. 

A larger sample is required to fully characterize the 
properties of these objects and probe the full spatial and 
velocity distribution of this population. Establishing the 
origin of these objects would provide important clues 
for understanding galaxy and cluster evolution. Objects 
which are remnant nuclei from a previously larger popu- 
lation of dE,N would be useful probes for understanding 
cluster destructive processes, and could help mitigate the 
missing satellite problem. Objects which prove to be star 
clusters would be valuable probes of galaxy merger and 
cluster merger and dynamical histories. 
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Table 1 
Keck/LRIS observations 



Mask 


Date 


a (J2000.0) 


8 (J2000.0) 


PA (deg) 


Seeing (arcsec) 


No!,j 


1 


2 Apr 2008 


13 00 40.69 


28 01 54.86 


1.5 


0.7 


35 


2 


3 Apr 2008 


13 00 17.82 


28 01 26.81 


1.5 


1.0 


32 


3° 


3 Apr 2008 


13 00 24.62 


27 56 26.11 


80.0 


1.0 


38 


4 


2 Apr 2008 


12 59 40.10 


27 59 06.28 


121.0 


0.8 


39 


5" 


28-29 Mar 2009 


13 00 43.95 


27 59 47.08 


85.0 


0.8-1.4 


43 


6 


30 Mar 2009 


12 59 53.57 


27 57 50.93 


-94.0 


1-1.3 


42 



^Exposure times for both red and blue chips were 8 X 1500 sec with the exception of Mask 3 for 
which we obtained only 6 X 1500 sec on the blue side. 



Observations were taken in thick and variable cloud cover. 
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Table 4 
Effective Radii 





GALFIT 


{Scrsic) 






ISHAPE 








ID 


Re (pc) 


n 


F814W 




profile 


index 


R. (pc) 





191006 


7.1 ± 2.3 


5.2 


22 


51 





82 


King 


30 


7.8 ± 2 


6 


0.68 


192636 


33.3 ± 3.8 


3.1 


22 


28 





86 


King 


30 


29.2 ± 6 


5 


0.93 


120985 


40.5 ± 3.2 


1.7 


22 


17 





90 


King 


30 


42.8 ± 12 


2 


0.88 


121666 


42.0 ± 13.3 


3.9 


21 


54 





99 


King 


100 


38.3 ± 3 


6 


0.88 


195526 


37.6 ± 14.9 


5.7 


21 


66 





91 


King 


100 


32.5 ± 5 


3 


0.99 


195614 


22.0 ± 2.5 


6.9 


22 


23 





84 


King 


100 


22.9 ± 5 


2 


0.89 


196790 


20.9 ± 2.4 


6.6 


21 


74 





86 


King 


100 


22.6 ± 6 


1 


0.87 


182204 


14.7 ±4.0 


5.4 


23 


17 





65 


King 


100 


8.6 ± 4 


6 


0.91 


242857 


23.7 ± 3.0 


4.7 


22 


30 





79 


King 


30 


24.8 ± 5 


2 


0.79 


160141 


7.0± 3.1 


5.1 


22 


66 





76 


King 


100 


6.0 ± 2 


9 


0.89 


161244 


4.8 ± 2.9 


3.7 


22 


60 


1 


00 












163341 


23.3 ± 2.7 


7.6 


22 


24 





98 


King 


100 


19.5 ± 6 


7 


0.87 


163400 


66.9 ± 16.4 


6.3 


23 


11 





74 


King 


100 


54.4 ± 7 


1 


0.80 


92415 


24.3 ± 3.8 


7.6 


22 


67 





91 


King 


100 


22.0 ± 6 


4 


0.87 


163575 


17.1 ± 4.1 


5.3 


22 


53 





74 


King 


100 


12.8 ± 6 


6 


0.90 


150880 


12.8 ± 5.9 


4.4 


22 


74 





95 


King 


100 


11.0 ± 2 


4 


0.77 


151072 


125.5 ± 18.8 


6.6 


21 


29 





98 


King 


100 


68.1 ± 6 


4 


0.89 


150000 


36.0 ± 8.3 


7.9 


22 


38 





92 


King 


100 


24.3 ± 2 


8 


0.72 


81669 


18.2 ± 6.0 


2.6 


22 


95 





95 


King 


30 


17.7 ± 5 


4 


0.97 


1041346 


27.7 ± 7.5 


4.4 


21 


89 





94 


King 


100 


28.3 ± 4 


5 


0.85 


1039188 


25.6 ± 4.2 


5.2 


21 


97 





89 


King 


100 


26.4 ± 6 





0.88 


1041508 


10.2 ± 3.9 


6.1 


22 


91 





74 


King 


100 


14.4 ± 4 


6 


0.86 


1043225 


20.9 ± 7.0 


4.2 


22 


92 





64 


King 


100 


13.0 ± 4 


2 


0.74 


1042830 


19.7 ±4.8 


4.1 


21 


87 





97 


King 


100 


15.2 ± 6 


8 


0.94 


2000005 


7.1 ± 3.3 


4.1 


23 


07 





85 


King 


100 


6.7 ± 2 


5 


0.98 


1044251 


34.7 ± 8.4 


4.5 


22 


15 





84 


King 


30 


32.6 ± 6 


7 


0.85 



1044847'' 



''This object is outside of our ACS coverage. Although this object has been observed in archival WFPC2 
images, the larger pixel scale makes size measurements for such small objects impossible. 

^Typical uncertainties in b/a were 0.07 and 0.05 for Galfit and Ishape, respectively. However, for such 
small objects, we don't consider the axial ratio measurements to be very reliable. 



Table 5 
Age and Metallicity 

Type N Hp Fe 5270 Mgb Age [Fe/H] [Mg/Fe] V - /" [Fc/Ii]^v - 1)^ 

AAA Gyr dex dcx mag dex 



UCDs- red 


6 


1.53 ± 0.23 


2.40 ± 0.32 


2.42 ± 0.29 




-0.62"J2 


-0.10°;;^ 


1.18 


-0.41 


UCDs- blue 


9 


2.41 ± 0.31 


0.99 ± 0.47 


1.33 ± 0.36 




0.95 


-1.40 


UCDs- bright 


5 


2.01 ± 0.34 


1.84 ± 0.45 


2.73 ± 0.40 


8.7^*^ 


-o.8o«:g 


23 






UCDs- red,N4874 


5 


1.73 ± 0.25 


2.22 ± 0.33 


2.60 ± 0.31 


12.6j-r 
4.00^:? 


1.18 


-0.40 


UCDs- strong H/3 


8 


2.99 ± 0.28 


1.09 ± 0.41 


1.98 ± 0.32 


-1.24—,- 


0.97 


-1.31 


UCDs- weak H/3 


8 


1.20 ± 0.28 


1.92 ± 0.39 


2.02 ± 0.39 


1.08 


-0.84 


UCDs- metal rich 


2 


1.63 ± 0.43 


3.08 ± 0.54 


2.52 ± 0.52 


13.2"- 
5 2 

7 4^-^ 
' •*2.0 


-0 

— 7^?^-^^ 
-0.660 16 


— 21 

(!?9 

— 07^.^ 

H4fi 
-0.24 ; 

-0.18ij;f^ 


1.17 


-0.47 


UCD/dE,N 


2 


2.48 ± 0.22 


0.76 ± 0.31 


1.03 ± 0.27 


1.07 


-0.87 


dE,N(brt) 
dE,N(fnt) 


4 
24 


2.26 ± 0.19 
2.42 ± 0.14 


1.85 ± 0.27 
1.84 ± 0.20 


1.28 ± 0.24 
1.04 ±0.17 


1.08 
1.06 


-0.82 
-0.91 


dEs 


10 


2.19 ± 0.29 


1.29 ± 0.39 


0.82 ± 0.33 


-1.16"J9 







Individual: 

150000 3.01 ± 0.62 1.57 ± 0.85 1.82 ± 0.75 2.1^;? -0.40",^" O-O^nio 

196790 1.95 ± 0.66 2.71 ± 0.91 4.18 ± 0.83 7.3^^^ -0.08S-|j 0.24^|^ 1.25 -0.12 

121666 1.84 ± 0.44 3.48 ± 0.57 4.27 ± 0.51 7.6;^^^ ^"^V^^ '^''^^RA^ ^'^'^ -0.47 

cE (91543) 2.08 ± 0.12 3.04 ± 0.15 4.22 ± 0.13 4.8f:| 0.08^-g9 0.20^;§^ 1.33 0.22 



Average values are used for combined spectra. 
^Estimated from the relation between (V — /) color and [Fe/H] provided by |Barmby et al.| | |20o'o| 

Note. — Where age and metallicity values arc absent, this is duo to measured line indices falling outside of model ranges. 
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Table 6 
Velocity Distribution 



Sample 


Number 




err 


(7^ 


err 






km s ^ 




km 




Full LSB 


51 


6970 


178 


1269 


126 


Full UCD 


27 


6887 


207 


1072 


146 


{V - I) > 1.05 


11 


6992 


379 


1260 


269 


(V - I) > 1.10 


8 


6907 


325 


922 


231 


(V - I) < 1.05 


14 


6707 


231 


861 


163 


a < 195 


12 


7296 


160 


558 


114 


a > 195 


15 


6559 


325 


1257 


230 


N4874 proximity 


9 


7257 


87 


254 


60 


N4889 proximity 


5 


6526 


148 


332 


105 
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Fig. 1. — Spatial distribution of the UCDs. Red boxes are the locations of observed ACS fields, black rectangles are the 6 LRIS masks. 
Green points represent our original candidate sample, brown points are the expanded candidate sample, and cyan points are candidates chosen 
strictly on the basis of color. Lar ge solid circles de note the location of confirmed UCDs and triangles as less certain UCDs. Open circles 
mark the location of compact dEs ( [Price et al.|200"9] |. X's are UCD candidates determined from redshifts to be stars (black) and background 
galaxies (red). 
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4000 4500 5000 5500 



wavelength (A) 

Fig. 2. — Spectra for 6 of our UCDs, smoothed 3 times. For each, we list the object ID as in Tables [2)3| along with the radial velocity in 
km s~^ derived from absorption lines. 



22 



Chiboucas et al. 




4000 4500 5000 5500 



wavelength (A) 

Fig. 3. — Spectra for another 6 UCDs, smoothed 3 times. Labels as in Figure [2] 
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Fig. 4. — Spectra for 7 confirmed UCDs, smoothed 3 times. Labels as in Figure|2] 
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Fig. 5. — Spectra for the second set of 8 confirmed UCDs, smoothed 3 times. Labels as in Figure[2] 
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Fig. 6. — ACS F814W— band thumbnail images, 17.5 arcsec across for all 19 confirmed UCDs. Circles that identify the UCDs are 2.6 arcsec 
in diameter. 
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Fig. 7. — ACS F814VF— band thumbnail images, 17.5 arcsec across for 8 UCDs in the second sample. The WFPC2 image containing 
1044847 comes from archival HST data (program G06283, PI J. Westphal). 
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Fig. 8. — Color-magnitude diagrams with UCD candidate selection criteria shown. Photometry comes from| Adami et al.| | [2006l l. Objects 
classified as stars and galaxies come from the full Adami et al. catalog. The original and extended UCD candidate samples are shown separately. 
Boxes bound the original color/magnitude cuts used to select UCD candidates within our ACS fields. dEs are Coma cluster memb ers observed 
with the same masks as the UCDs, while open circles represent cE member galaxies with redshifts measured from Hectospec data l |Price et al.| 
|2009| Marzke et al. in prep.). 
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Fig. 9. — SExtractor FWHM vs. magnitude (left) and SExtractor classification vs. FWHM (right). FWHM is in pixels, measured from 
our ACS data. Small pentagons are objects from our first observing run with spectra having too low S/N to measure redshifts. 
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Fig. 10. — Colors for all objects in our spectroscopic sample. BVRI colors come fromlAdami et al.J p006| l, F475VK - F814W {g - I) from 
our ACS imaging. Magnitudes are extinction corrected using the dust maps of Schlegel et al. (1998). Light circles are objects with a faint 
but obvious LSB envelope around a nucleus. Other symbols as in Figure [s] Representative errors for the UCDs are provided for the ACS 
photometry {g — I vs. /) and Adami et al. photometry {B — V vs. R). Photometry docs not exist for all objects in all bands, so a few UCDs 
are missing from various panels in this figure. 
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Fig. 11. — Left: central surface brightness vs. ij— band total magnitude for all objects in our ACS survey region. Photometry comes from 
[Adami et al.| (2006). Ri ght: mean effec tive surface brightness vs. ^81414^— band total magnitude for all objects in our ACS survey region. 
Fhotometry comes from [Hammer et al.| l|2010[| . 



30 



Chiboucas et al. 




20 40 60 



GALFIT (pc) 

Fig. 12. — Comparison of Galfit and Ishape measured Re. One object, 151072, has a much larger measured size and is not shown here but 
we note that the Galfit measurement for this object deviates significantly from that of Ishape. 
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Fig. 13.— Residuals from best Galfit (middle) and Ishape (right) fits for 4 UCDs. Prom top to bottom: 150000, 121666, 195526, and 
151072. Panels are 2.5 arcsec ax;ross. Strong residuals are apparent and may indicate that two-component fits are warranted in several of 
these cases. 




Fig. 14.— Residuals from best Galfit (middle) and Ishape (right) fits for another 4 UCDs. Prom top to bottom: 1041508, 182204,195614, 
and 163400. Panels are 2.5 arcsec across. 
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Fig. 15. — Measured Fe 5270A vs H/3 line strengths. [Fe/H]-age grid models are from Schiavon 
Plotted are Fe 5270 and H/3 measurements for composite spectra of similar object types 



1 2007 1 for a/Fe = (blue) and 0.3 (cyan). 



Object types are labeled on the plot. We include 
two points for dE,N; the larger solid circle includes 4 objects with bright, prominent nuclei, the smaller one includes 24 dE,N with small, faint 

2009) for the same object labeled as cE below it. Star-like symbols refer to UCDs with 
different combinations of stacked spectra:" the 5 brightest UCDs, 6 red (V - / > 1.05), 9 blue {V - I < 1.05), 5 red UCDs around NGC 4874, 

The UCD/dE,N includes object 151072 from this paper and 242439 from |Chiboucas| 
n extended envelope surrounding them. We also plot 3 individual bright UCDs (small 



nuclei. The open black circle comes from Price et al 
tit combinations of stacked spectra: the 5 brig 
8 with weak H/3{< 2.5), and 8 with strong H/3(> 2.5) 
Both are compact sources with 
not include the large error bars 



|et al. l |2010| . Both are compact sources with a hint of an i 
stars j but do i 
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Fig. 16. — Histograms of UCD radial velocities. Left: UCD radial velocity distribution compared to that of dEs with membership 
spectroscopically determined from the same LRIS MOS observations. The location in velocity space of several prominent giant galaxies are 
labeled. Middle: UCDs separated into blue and red populations by V — / color. Right: UCDs split by RA. 




Distance from IC 4051 (arcmin) Distance from NGC 4874 (arcmin) 

Fig. 17. — UCD radial velocity as a function of projected distance from prominent Coma cluster giants (asterisks) including the cD galaxies 
NGC 4874 and NGC 4889, and IC 4051 (a galaxy with one of the highest Coma cluster globular cluster specific frequencies). One arcmin 
corresponds to 29 kpc. Red symbols have (V — I) > 1.05, while blue have {V — I) < 1.05. Where V — I colors do not exist, we use colors 
in other bands to infer a rough {V — I) color. Larger symbols represent brighter magnitudes. Open circles are cEs, triangles are objects 
with insecure redshifts. Objects which may be associated with individual galaxies are encircled. The bottom right panel displays all possible 
associations of UCDs with giant galaxies. 
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Fig. 18. — Left: Cumulative distribution of the separation of each UCD from the nearest of one of 3 giants in the core region: NGC 4874, 
NGC 4889, and IC 4051. 25,000 points randomly distributed within the LRIS footprint are generated to produce the expectation for a uniform 
distribution. Right: Cumulative distribution of the velocity difference between each UCD and the closest of the 3 giants. The curve labeled 
'random' displays the same for 25,000 points with randomly sampled velocities assuming a Gaussian distribution with mean 6925 km/s and 
dispersion 1000 km/s. 
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Fig. 19. — Maps of confirmed and candidate UCDs (circles and squares, respectively), cEs (open circles), and giant galaxies (crosses) along 
with Top: candidate globular clusters (points ( |Peng et al.|20lH ) Bottom: confirmed and candidate dE (triangles) and dE,N (small circles) 
Coma members (Trentham et al. in prep). 
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Fig. 20. — Cumulative distribution as a function of distance from NGC 4889 out to tlie extent of tfie central ACS footprint. Top: projected 
distance from NGC 4889, Bottom: distance in declination only. The expectation for a uniform distribution (gray solid line) is determined by 
randomly distributing 10,000 points over the observed central ACS fields. 
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Fig. 21. — Color distribution of Coma cluster candidate globular clusters (points: full ACS survey region, small circles: visit 19 only) and 
confirmed UCDs (large circles: all UCDs, encircled: only those from visit 19). The dashed histogram displays binned _F814Vl^ > 24.7 visit 19 
globular cluster counts, the solid histogram includes only F814W < 24.7 visit 19 counts. We show the best double Gaussian fit to this set of 
brighter counts. The two peaks are at g — / = 0.93 and 1.17. 
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Fig. 22. — Magnitude-size relation for compact objects. Data for UCDs come from this work and 'Mieske ct al.' ("2002' '2007a' '2008bt 



Evstigneeva et al. (2008); Hau ct al. (20 09[|. Stru ctural parameters for other object types come from: cEs (Price et al. 2009 Kvstignccva ct'aL 

gan ct al." 



|Evs ^ 

[2008 Kent 1987), UG'i'Os (Ha§cgan ct al."2005), dE nuclei/nuclear star clusters ('Geha ct al. 2002 Gcorgiev ct al. 2009 Cote ct al. 201 
othe r intermediate compact types (Mieske et al. 2008b, 2007a) , globular clusters (Fusi Pccci ct al. 1994 Grillmair et al. 1996, Barmby et"al 
|2002| |McLaughlin fc van der Maren[20J5| jPa Gosta et al.||2[TO| , and Goma clusteF'd Es (Ghibo ucas et al. 2010 Hoyos ct al.|201Tj i: DasEe; 
lines represent lines ot constant surface brightness. We use the relation from |Peng et al. ( 2006 1 to convert jGotc et al._ ( ^20d6} g'-band data 
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Fig. 23. — Size vs color and distance from NGC 4874 vs size for confirmed UCDs. Best fit linear relations, excluding two outliers, is shown 
as the dashed line. 



